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SUMMARY 

A model to predict charging effects in the Jovian magnetosphere was developed for 
the preliminary design of a Jupiter probe. Charging calculations using this model are 
presented and discussed. At distances from Jupiter of less than 14. 6 Jovian radii, dif- 
ferential charging of interior structures relative to adjacent exterior s\xrfaces may be as 
much as 15. 4 kilovolts. Kilovolt differentials are also predicted between the sunlit and 
dark exterior surfaces of the probe and in the region of the Jovian magnetosphere where 
corotation of the low-energy plasma affects ion collection. Design modifications to re- 
duce these differentials are recommended. 


INTRODUCTION 

Electrostatic charging of spacecraft surfaces to kilovolt potentials occurs when 
clouds of energetic particles are injected into the midnight region of the magnetosphere 
near synchronous altitude during geomagnetic substorm activity. The first evidence of 
spacecraft charging to kilovolt potentials was reported by DeForest (ref. 1) in 1972, 
based on data from the University of California at San Diego (UCSD) plasma detector on 
the ATS -5 spacecraft. Since that time, a number of workers have studied aspects of the 
spacecraft charging phenomena. Grard (ref. 2) and Rosen (ref. 3) provide extensive in- 
formation on the subject. Existing knowledge of spacecraft charging was used in this 
study to develop a theory that would predict charging for a Jupiter probe in the Jovian 
magnetosphere. 

Any svirface immersed in a charged -particle environment will come into equilibrium 
with it by acquiring a potential such that in steady state the net cvirrent to the surface is 
zero. Since for most plasmas the electron flux is much larger than the ion flux, the sur- 
face acquires a negative potential. Most materials respond to charged -particle bombard- 
ment by backscattering or reflectii^ electrons and by emitting secondary electrons. 

These emitted particles effectively reduce the electron current to the siirface and thus 



reduce the negative potential. In addition, illuminated sxirfaces emit photoelectrons, 
which contribute an additional outgoing electron current. 

Thus it is not sxirprising that spacecraft can become charged by their environments. 
What was surprising in the geosynchronous case was that the particles injected into the 
m^netosphere by geomagnetic substorms have kilovolt energies and thus can charge 
spacecraft surfaces to kilovolt potentials. 

The Jovian magnetosphere has large fluxes of energetic particles. Scarf and 
Fredericks (ref. 4) suggest that these particles may interact with spacecraft, possibly 
charging their surfaces to large potentials. In addition, because many of these particles 
have sufficient energy to penetrate into the interior of a spacecraft, a spacecraft's in- 
terior structures may become charged relative to its exterior surfaces if they are 
electrically isolated from each other. This "inside out" charging is possible for the 
Jupiter probe since the proposed exterior surface is entirely insulating and its metal 
structure has no electrical access to the exterior. 

Because of concern about possible charging effects on the Jupiter probe's perform- 
ance, an analytical study was performed. The specific objectives of the study were to 
investigate the possible charging of both the exterior and interior surfaces of the probe 
during its descent through the Jovian m^netosphere and to make design recommenda- 
tions to minimize identified chargii^ problems. This report describes this chargii^ 
study and presents and discusses the results. 

Jon C. C^lebay of the Lewis Research Center wrote, debugged, and ran the com- 
puter program used in this study. 


JOVIAN CHARGED-PARTICLE ENVIRONMENT 

A model of the charged -particle environment in the Jovian magnetosphere, based on 
data from the Pioneer 10 and 11 spacecraft (ref. 5), has been developed by the Jet 
Propulsion Laboratory. This model gives fluxes of high-energy particles in the Jovian 
magnetosphere and the characteristic temperatures of the magnetospheric plasma com- 
ponents: ~4 eV for electrons, ~100 eV for protons inside about 6 Jovian radii, and 
~400 eV for protons from 6 radii to the magnetopause (~60 radii). These low-energy 
plasma components are assvimed to be Maxwellian and are hereinafter called thermal. 
Because they are trapped by the planet's magnetic field, they corotate with Jupiter. 

The fluxes of high-energy particles are considered to be isotropic in this model. 
They are not expected to be Maxwellian and do not corotate with Jupiter. In this study, 
"high energy" particles are electrons with energies from 0. 06 to 100 MeV and ions with 
energies from 1 to 100 MeV. The limits on these energy ranges are the lowest and high- 
est energies for which information was available. Figure 1 shows the fluxes of thermal 
electrons (e > 0) and protons (e > 0) and of high-energy electrons (e > 0. 06 MeV) and 

2 



protons (Cp > 1 MeV) in this study as a function of time to Jovian encounter along a pos- 
sible mission trajectory. The long ticks on the abscissa indicate times for which probe 
potentials were calculated. 

In the present study of Jupiter probe chargii^, the high-energy and thermal (low 
energy) cvmrents are modeled separately. Their values are then combined to achieve 
the required ciirrent balance because the dependence of the particle currents on the probe 
surface voltages is expected to be different for the low- and high-energy particles. 


JUPITER PROBE MODEL 

The Jupiter probe configuration studied here has a blunted, 60° half -angle, conical 
forebody and a hemispherical afterbody. The forebody is a carbon-phenolic, ablative 
heat shield that protects the aluminum interior structure and instrumentation. The 
afterbody is a fiberglass -phenolic, honeycomb hemisphere (45-cm spherical radius) that 
protects the antenna. The outer surfaces of both the forebody and the afterbody are 
covered with multilayer -aluminized-Mylar thermal blankets. For the charging model, 
the probe is considered a sphere of 45 -centimeter radius. Figure 2 shows the probe 
configuration and the spherical model. 

The probe is separated from the orbiter 57 days before Jupiter encoimter (at about 
500 Jovian radii) and coasts toward the planet on a near -zero -inclination trajectory at a 
nominal spin rate of 5 rpm about its axis of symmetry. During the coast, the afterbody 
is oriented toward the Earth. As the probe enters the sensible atmosphere (altitude of 
~450 km), it tips about 90° so that the forebody heat shield points toward Jupiter and the 
afterbody toward the orbiter. Nominal velocity before entry is about 60 km/sec. The 
probe and the mission are described in detail in references 6 and 7. 

Since the coast phase of the mission is of primary interest for the charging study, 
the relative geometry assumed in the model approximates that of the coast phase. This 
geometry is shown in figure 3. Because the Sun is assumed to be in the same direction 
as the Earth relative to the probe, the afterbody is sunlit and the forebody, dark. In ad- 
dition, the probe’s velocity vector is assumed to be perpendicular to the radial line con- 
necting the center of Jupiter to the center of the probe. In this way the effects of the 
relative motion of the probe and the thermal plasma - which, as noted in the previous 
section, corotates with the planet - on particle current collection can be more easily 
calculated. 

As a quick review of the probe's structure, the exterior surface is completely cov- 
ered by insulating thermal blankets, the entire afterbody is an insulating structure, and 
the conducting aluminum structure is in the forebody. If the almninized layers of the 
thermal blankets are not groimded to the aluminum structure, the only mechanisms for 
charge transfer from the forebody to the afterbody are leakage along the insulator 
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surfaces or leakage through the sheath. The present model treats the forebody and 
afterbody as independent, an assumption that affects the predicted voltage differences be- 
tween the forebody and afterbody surfaces as discussed in the results section. 

The model currents are defined in figure 4. The equilibrium potentials for the four 
surfaces of interest are for the exterior surface of the forebody, for the alum- 
inum structure that forms the interior surface of the forebody, cp for the exterior sur- 
face of the afterbody, and a for the interior surface of the afterbody. The currents 
^HE f %E a caused by high-energy particles impinging on the forebody and 
afterbody exterior surfaces, respectively. These particles include electrons and pro- 
tons that are incident on the surfaces but do not have sufficient energy to completely pen- 
etrate the appropriate insulating layers, as well as the backscattered and secondary 
electrons generated by them. The net cvirrents I-pjj j and ^ are caused by ther- 
mal plasma particles impinging on the forebody and ^terbody exterior siirfaces, respec- 
tively. These particles include both primary and secondary (emitted) particles. The 
photocurrents IpjjQ j and IpjjQ ^ are emitted by the respective surfaces. The leak- 
age currents through the bulk of the insulators are Ij^ £ and The current 

Ihep a caused by high-energy protons and electrons that have penetrated to the alu- 
minum structure through either the forebody or afterbody insulators and are stopped by 
the aluminum structure. In addition, those particles coming through the afterbody in- 
sulators that have insufficient energy to penetrate completely through the probe and es- 
cape contribute to this current. Thus, Ij^gp ^ is the particle current directly incident 
on the aluminxun structure. The current Ijjgp £ is caused by high-energy particles 
that have penetrated both the forebody insulators and the aluminum structure but cannot 
penetrate completely through the probe and escape. Thus, Ijjgp £ is the particle cur- 
rent directly incident on the interior surface of the afterbody. ’ 

In summary, the charging environment for the exterior surfaces consists of cur- 
rents caused by high-energy particles, thermal plasma particles, and photons. The 
charging environment for the interior simfaces consists of currents caused by high- 
energy particles. The interior and exterior surfaces are connected by leakage paths 
through the insulating layers. 


PROBE CURRENT BALANCES 

Thermal Plasma and Photoelectron Currents 

The thermal plasma currents to the probe include contributions from primary elec- 
trons and ions and from secondary electrons due to electron impact. That is, they are 
of the form 
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( 1 ) 


where cp^ is the appropriate surface voltage and 1^, 1^, and Ig are the primary elec- 
tron, primary ion, and secondary electron cxirrents, respectively. The ion cvirrents 
also depend on whether the probe is at rest or in motion with respect to the corotating 
plasma: At distances nearer the planet than about 23 Jovian radii, the probe is effec- 
tively at rest; further from the planet, the probe is in motion. That is, at 23 Jovian 
radii or less, the ion thermal velocity is greater than the relative velocity of the probe 
and the corotating plasma; at 23 Jovian radii or more, the relative velocity exceeds the 
ion thermal velocity. 

Table I summarizes the functional forms for the thermal and photoelectron currents 
to the probe, for the geometry of figure 3. Figure 5 shows the thermal and photoelec- 
tron cxirrents to the probe. These currents and their functional forms are discussed in 
more detail in appendix A. 


High -Energy -Particle Currents 

Fluxes of high-energy particles can charge both the exterior and interior surfaces 
of the probe since some particles will penetrate into the interior. The currents caused 
by high-energy particles impinging on the probe should not be affected by the probe’s sur- 
face voltages unless these vpltages become comparable to the particle energies. These 
currents are therefore considered to be independent of until |e9>g| becomes approx- 
imately 0, 06 MeV. 

The information on high-energy particle fl\ixes used in the present study is in terms 
of integral fluxes of electrons and ions along the probe's trajectory. Currents caused by 
particles with energies greater than can be foxmd from expressions of the form 

I = iAer(€^) (2a) 

a 

where A is the surface area of interest, e is the electronic charge, and F(cj) is the 
integral flux at Cj (i. e. , the flux of particles with energies equal to or less than ej^). 
Currents caused by particles with energies between €j and Cg (^2 ^ found 

from expressions of the form 


I = iAe[r(€i)-r(c2)] 


(2b) 
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The factor of 1/2 in each of these equations accounts for the fact that the model used to 
calculate the fluxes assumes an isotropic particle distribution: Only half the particles 
will have velocities directed toward the probe. 

Figure 6 depicts the behavior of high-energy particles impinging on the probe. Par- 
ticles impinging from the forward direction but unable to completely penetrate the fore- 
body insulating layers (e < e^) and those imputing from the aft direction but \mable to 
completely penetrate the afterbody instilating layers (e < e„) are considered to contribute 

a 

to the currents to the forebody and afterbody exterior simfaces, respectively. 

Particles impinging from the forward direction with energy or greater but un- 
able to completely penetrate the aluminum structure and enter the afterbody region 
(cj ^ e < and particles impinging from the aft direction with energy or greater 

but unable to completely penetrate the probe and escape (e„ s e < e ) are considered to 
constitute the currents to the internal aluminum structure. All particles approaching 
from the aft direction that would be stopped either in the aluminum structure itself or in 
the forebody insulators are considered to constitute a current to the aluminum structure, 
which is effectively the interior surface of the forebody. Those particles approaching 
from the forward direction with energies ^ or greater but unable to also penetrate 
the afterbody insulators and escape (e^ ^.1 “ ^ ^e^ considered to contribute to the 

current to the inside of the afterbody. Finally, particles from either direction with suf- 
ficient energy to penetrate completely through the probe and thus escape (e s are not 
considered sources of charging current to any part of the probe. 

These assxxmptions about the high-energy particle currents ignore the effects of 
charge deposition in the insulator bulk. Although consideration of such effects is beyond 
the scope of this study, it is warranted in further work. Each of these energies (cj, 

€f etc. ) is different for electrons and ions since ions have much smaller penetrating 
ranges in materials than electrons of the same energy. This is accounted for in the 
model. Expressions used for the various high-energy currents are summarized in 
table n. Their functional forms are discussed in more detail in appendix B. 


Current Balances 

The cxirrent balance is solved to determine the potentials as follows: Balancing the 
forebody current requires that the net currents to both the heat -shield exterior and the 
aluminum structure be zero. Thus, 

® ^ ■ VHO,f^'*0f) " ^L,f^^f’ ^Al^ 


and 
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(4) 


° “ %EP, aW 

The leakage current through the shield depends on the difference between <pj and <p^ 
so that 


, _ <Pf- <Pm 

Rj 


(5) 


where Rj is the resistance of the forebody insulators. Eliminating ^ from equa- 
tions (3) and (4) yields ’ 


° " %E , ^H, f (*i^f ^ ■ ^PHO, f (^f ^ ^HEP , a^^i^f ^ 

which can be solved for <p^. ^ Then is calculated from equations (4) and (5), 

Similarly, balancing the afterbody cttrrent requires that the net currents to both the 
interior and exterior surfaces of the afterbody be zero. The equations, analogous to 
equations (3) to (5), are 

° = ^HE, a^*i^a) ^TH, a^^a^ " %HO, a^*i^a^ ' ^L, a^'^a» ^a.r^ 

0 = lHEP,f(<Pa^ + ^L, 



<P^ - (P, 


ar 


R. 


(9) 


where R„ is the resistance of the afterbody insulators. After is found from equa- 

A -4 ^ A 

tions (7) and (8), ^ <p„_ can be determined from equations (8) and (9). 

a£ 


Specific Probe Parameters 

Solving the current -balance equations requires assigning numerical values to a nxim- 
ber of parameters that depend on the materials and construction of the probe: stirface 
areas, resistances, secondary emission and backscattering coefficients, penetration en- 
ergies for the high-energy particles, etc. Table IE shows each parameter, its symbol, 
the value used in this model, and the basis on which this value was chosen. Values for 

^The currents are expressed as functions of the potentials in tables I and n. 


7 



the backscattering and secondary emission coefficients were, as indicated, estimated 
from ranges or estimates given in references. Values for the resistance and the pen- 
etration energies were calculated from several pieces of information (appendix C). 


RESULTS AND DISCUSSION 

The ciirrent -balance equations were used to calculate the equilibrium potentials of 
the forebody exterior surface cp^, the aluminiim structure cp^, the afterbody exterior 
surface cp^, and the afterbody interior siirface (p^ for the positions along the probe’s 
trajectory indicated in figure 1. The transition from "in motion" to "at rest" (at 23 
Jovian radii) takes place about 24. 5 hours before Jovian encounter. (See appendix A for 
further discussion of this transition. ) 


Forebody Potentials 

Potentials predicted for the forebody exterior surface <p^ and the aluminum interior 
structure <p^. are shown as functions of position and time to Jovian encounter in fig- 
ure 7. Significant charging of the forebody exterior surface with respect to the plasma 
occurs from 7 to 14 hours and more than 24 hours before encounter. These times corre- 
spond to distances from the planet of 9 to 15 Jovian radii and more than 23 Jovian radii, 
respectively. Because in the latter region the probe is in motion with respect to the co- 
rotating plasma, the ion component of the thermal plasma current is cut off in this model. 
Because the forebody exterior stirface is in the dark, the only source of ion current to it 
is the high-energy ions. The flux of these ions is very low (fig. 1). There is charging 
from 7 to 14 hours before encounter probably because the thermal plasma is simply in- 
sufficient to balance the effects of the high-energy electrons. The most negative poten- 
tial (p^ predicted here is -18 kilovolts at 12 hours, or 13 Jovian radii. These predicted 
potentials are overly pessimistic in regions where large differentials are expected be- 
tween the forebody and afterbody exterior surfaces, as discussed on page 10. 

The almninum structure is also predicted to charge in two regions: More than 24 
hours before Jovian encounter, <pj^ is about equal to (p^. Less than 14 hotirs before en- 
counter, (pp^ is more negative than <p^, <pp^ reaches its most negative potential 
(-21 kV) at 11 and 12 hours. In addition to this "absolute" chargii^ (with respect to the 
plasma), there is significant differential charging between the aluminum structure and 
the forebody exterior surface less than 14 hours before encounter. This aspect is dis- 
cussed in the section Differential Charging. 
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Afterbody Potentials 


Potentials predicted for the interior cp^ and exterior afterbody surfaces are 
shown as functions of position and time to Jovian encounter in figure 8. The afterbody 
exterior surface experiences almost no chargii^ and is generally a few volts positive. 

No significant changes occur in the transition from motion to rest. The photocurrent ap- 
parently maintains the potential very near plasma groimd. The afterbody interior shows 
some charging less than 9 hours before Jovian encounter, or less than 11 Jovian radii, 
but the charging is much less severe than for the forebody exterior surface. Its maxi- 
mum negative potential is -195 volts, which is not very significant from a charging 
standpoint. 


Differential Charging 

Figures 7 and 8 indicate two possible problems in differential charging: the differ- 
ential between the forebody exterior surface and the aluminxim interior structure and 
that between the forebody and afterbody exterior surfaces. There is also a differential 
between the aluminvim interior structiire and the afterbody interior surface, but solving 
the first two problems should automatically solve it as well. 

Figure 9 shows the differential potentials - <p^ and (p^ - (p^ as functions of 
time to enco\mter. It is notable that <pp^ is always more negative than or equal to (p^ 
and that <p^ is more negative than cp^. 

Differential across forebody insiilators. - There is a kilovolt potential across the 
forebody insulators (fig. 9) less than 13 hours (14. 6 Jovian radii) before encounter, with 
a maximum of 15. 4 kilovolts at 9 hovurs (11 radii). The probe forebody is composed of 
the Mylar thermal blankets, the carbon -phenolic heat shield, and a fiberglass structure 
(fig. 13). In obtaining the overall resistance of the forebody insulator, these components 
were treated as resistors in series (appendix C). If these resistors behave as a voltage 
divider, the 15.4 kilovolts across the insulators yields 8. 7 kilovolts across the Mylar, 

— fi 

6. 7 kilovolts across the fiberglass, and ~5xl0" volts across the heat shield. These 
voltages indicate stresses of ~3xl0 volts per meter (~750 V/mil) on the Mylar and 
~3xl0^ volts per meter (~7. 5 V/mil) on the fiberglass. Since the dielectric strength of 
Mylar is several thousand volts per mil, this stress should not cause bulk breakdowns of 
the Mylar. However, it may present a problem at edges or corners, where field inten- 
sification is expected, and which are thus the most probable regions for arc breakdowns 
to occur. The present model cannot treat edges and corners. Assessing the fields in 
such regions requires a more detailed model of the sxirfaces. 

These voltage stresses probably would present a serious problem. However, since 
there is differential charging, care should be taken to consider its possible effects. For 
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example, simply tying the aluminized layers of the Mylar to spacecraft ground, in the 
standard fashion, may cause all the differential potential to appear across the external 
Mylar film and thus approach and perhaps exceed Mylar's dielectric strength. Two sim- 
ple ways to minimize the differential potentials in the forebody insulators might be to 
metalize and groimd the exterior of the Mylar or remove the Mylar entirely, and then in 
each case to install a grounded conductive layer between the fiberglass and the heat 
shield or otherwise provide a conductive path around the fiberglass to the aluminum 
structure. ("Ground" here refers to attachment to the aluminum structure. ) 

Differential between forebody and afterbody exteriors. - The potential difference 
between the forebody and afterbody exterior sxirfaces is shown as the dashed line in fig- 
ure 9. Differentials to 18 kilovolts are predicted from 7 to 14 hours before encoimter 
and differentials to 300 kilovolts, more than 24 hours before encounter. These predic- 
tions are based on spherical probe theory and the assumption of ion shadowing more than 
24 hours before encounter. The negative hemisphere of a strongly differentially charged 
sphere can collect more ion current than one hemisphere of a sphere uniformly charged 
to the same large negative potential. Thus the potentials predicted for the forebody are 
more negative than if a more accurate asymmetrical model were used. Estimates of the 
effects of asymmetry indicate that kilovolt differentials can still be expected. From 7 to 
14 hours before encounter, the maximum differential is estimated to be at least 9 kilo- 
volts; more than 24 hours before encounter, a few kilovolts (rather than 300 kV) is esti- 
mated. Clearly, in reality the differential potential will not appear as a step change at 
the forebody -afterbody bovindary, since effects such as surface leakage and charge trans- 
port in the photoelectron sheath will assure a smooth transition. If the differential is 
distributed uniformly aloi^ the surfaces from forward to aft, the semicircumference n-rp 
of 1. 41 meters indicates fields of ~1. 3x10^ volts per meter along the surface for the 18- 
kilovolt total differential and ~10 volts per meter for the 3x10 -volt differential. If the 
differential is distributed other than linearly, there will be areas in which the stresses 
will be even larger. How the photoelectron sheath will behave in transferring charge 
from forward to aft is beyond the scope of this study, but it bears further investigation as 
the photocurrent is the critical factor in maintaining the afterbody surface potential near 
zero. 

Making the forebody surface relatively conducting would reduce the forebody - 
afterbody differential by allowing charge from the forebody surface to migrate freely. 

The forebody shield should also be electrically connected to siinlit areas so that photo - 
electron emission can remove electrons and reduce the negative charge. A photoelec- 
tron sheath could extend far enovigh to perform this function to some degree. But the 
best solution is to make the afterbody exterior surface conductive and connect it elec- 
trically to the forebody siuface. If this is not possible, a conducting annular structure 
electrically connected to the forebody siirface and exposed to sunlight should be 
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considered. However, any such structure will collect additional charged particles as 
well as additional photons. 


CONCLUSIONS AND RECOMMENDATIONS 

The Jupiter probe, as presently conceived, may have problems related to differ- 
ential charging since kilovolt potential differences develop between interior and exterior 
structures and between sunlit and dark surfaces. Although the voltage stresses across 
the bulk of the forebody insulators are predicted to be well below the bulk breakdown 
voltages of the insulators, edges and corners may present a problem. These voltage 
gradients can be minimized by incorporating relatively simple design changes. The dif- 
ferential between light and dark surfaces has less clear effects since their severity de- 
pends on the actual voltage distributions and material strengths. Again, it seems best 
to minimize such effects, and a simple design modification appears feasible. 

The following modifications are recommended to minimize differential charging: 

1. The forebody exterior surface should be made more electrically conductive by 
either removing the Mylar or metalizing the exterior surface. If the Mylar is metalized, 
all metalized layers of the thermal blankets, including the surface layer, should be 
grounded to the aluminum structure. 

2. An electrically conductive path should be created around the fiberglass, connect- 
ing the inside of the carbon -phenolic heat shield to the alximimun struchire. The after- 
body surface should be made electrically conductive; or a conducting structure, perhaps 
a metal band or rii^, should be added to give the shaded (now conductive) forebody sur- 
face electrical access to photoelectrons in order to reduce the light -dark differential 
voltages. 

With these modifications, particularly with the electrically conductive afterbody siirface, 
differential chargii^ should be reduced to a minimum. This is true even if the environ- 
ment is different from the one used in the model. The modifications do not, however, 
guarantee that the probe as a whole will remain near plasma groimd. If the probe must 
be maintained near plasma grovind (e.g. , for preentry scientific experiments), some type 
of potential -control device such as an electron emitter may be needed. 

This study also indicates some areas in which additional analysis would be useful: 
voltage as 3 numetry effects on particle collection; charge transport between the sunlit and 
dark svirfaces, both through the sheath and along surfaces; electric fields in the after- 
body cavity; and charge trapping in the insulator bulk. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 30, 1978, 

506-23. 
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APPENDIX A 


CURRENTS DUE TO THERMAL PLASMA FLUXES AND PHOTOELECTRONS 

For pxirposes of calculating the thermal currents to a Jupiter probe, the probe is 
assumed to be smaller than a Debye length so that orbit -limited theory applies. The 
thermal plasmas are assumed to be Maxwellian and to corotate with the Jovian magneto- 
sphere, an assumption that affects the ion collection in regions where the probe velocity 
relative to the plasma is greater than the ion random thermal velocity. Electron col- 
lection is not affected. 


Primary Particle Currents 

Electrons. - For a spherical collection geometry, currents due to collection of in- 
cident primary electrons are given by 


I.. = exp 
e ■’eo ^ 



for (Pg 



eo 



for cp„>0 

O 


(Ala) 


(Alb) 


where A is the surface area; cp^, the sxirface potential; e, the electronic charge; and 
kTg, the electron temperature. Here, is the random thermal cimrent density, which 
is given by 


'eo 


N e 
e 




(A2) 


where N^ is the electron number density; M^, the electron mass; and F^^, the inte- 
gral flux of electrons. 

Ions. - Since the thermal plasma is corotating with Jupiter, its corotation velocity 
is given by 


V - WjR = 12. 5 Rj km/sec 


(A3) 



where Wj is the angular frequency of Jupiter's rotation; R, the distance from Jupiter's 
axis of rotation; and Rj, the distance from the axis of Jupiter in Jovian radii. This 
velocity and the random thermal velocities of ions are plotted in figure 10 as functions of 
distance from Jupiter, The corotation velocity v exceeds the thermal velocity of ions 
at distances from the planet greater than about 23 Jovian radii. 

In regions where the probe is essentially at rest with respect to the plasma (i. e, , at 
less than 23 Jovian radii), the ion currents are foimd from 


Ii 2 = JioA exp 



for <pg -0 


(A4a) 


'l2 = 'lo^ 



for <p < 0 


(A4b) 


where kT- is the ion temperature and the ion random cxirrent density defined by 


■>10 


N.e 




(A5) 


where Nj is the ion number density; the ion mass; and r.^^, the integral flux of 
ions. 

In regions where the probe is in motion with respect to the plasma (i. e. , at 23 
Jovian radii or more), ion currents to the forebody and afterbody surfaces will be dif- 
ferent. This situation is illustrated in figure 11, Whipple (ref. 8) gives the ion current 
to the forebody surface as 



for > 0 


(A6a) 


/ 



< X 


D 


(A6b) 
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When cp is sufficiently negative that condition (A 6b) is not satisfied, 


I 


il 


N^ea-TTr 


2 

P 


V \<ii/\kT. „2 2 

\ l/\ 1 O'. / 



(A6c) 


In these equations, Vp is the probe radius; Xj^, the Debye length; the ion thermal 
velocity; and V, the relative velocity of the plasma and probe. For the geometry de- 
picted in figure 3, these forebody surface ion currents will impinge on the afterbody of 
the probe, and V will be given by 


V = V - 59 km/sec (A7) 

where v is defined in equation (A3) and 59 km /sec is the probe’s velocity. For an elec- 
tron plasma component having kT„ = 4 eV and N = 1 per cm*^ and for H'*’ ions, condi- 

C 6 r* 

tion (A6b) will be satisfied unless (p becomes more negative than about 4x10^ volts. 

O 

Since the afterbody is also svinlit (fig. 3), such a high negative potential is very unlikely. 
Therefore, equation (A6a) has been used in the present model. 

When the probe is in motion with respect to the plasma, there will be a relatively 
large void of ions "behind" the probe so that the side that is not being bombarded by the 
corotating ions (in this case the forebody shield) will effectively see no ion current. 


Emitted Currents 

Particle currents may be emitted by photoemission, reflection and backscattering of 
primary electrons, reflection of primary ions, and emission of secondary electrons due 
to ion and electron impact. From Whipple’s analysis (ref, 8), for 4-eV electrons and 
100- to 400 -eV ions the effects of reflection (for both electrons and ions) and of 
secondary -electron emission due to ion impact should be small (less than about 5 percent 
of primary ciirrents). No backscattering of electrons is expected for primaries with en- 
ergies less than 100 eV. Thus the emitted currents of importance here are photocurrent 
and current due to secondary electrons released by electron impact. 

Photoelectrons. - From the estimated density (1. 8x10 A/cm ) of a typical photo - 
current from a siirface at Jupiter, the sunlit half of the probe (the afterbody) will emit 
a photocurrent of 


IpHO = 1- 8X10"^° A^ for 0 


(A8a) 
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(A8b) 


IpHo = 1. 8X10-10 ^ j ° 

where A and <p are, respectively, the area and potential of the afterbody s\irface; 

A Ci 

and is the temperature of the emitted photoelectrons. In the present analysis 

is assumed to be 2 eV. The exponential term in equation (A8b) accounts for the 
reattraction of emitted photoelectrons to a positive surface. 

Secondary electrons released by electron impact. - The yield of secondary electrons 
released by electron impact is a function of the energy of the primary electrons and of 
the material properties. A functional form for this yield, attributed to Sternglass, has 
been given by Grard, et al. (ref. 9). The material properties needed are the maximum 
yield 6^^^ and the primary impact energy at which the maximum yield is attained. 
Willis and Skinner (ref. 10) give 6^^ =4.8 and = 175 eV for Mylar, which con- 
stitutes the outer sirrface of the probe. The yield c\irve and the fimctional form for the 
yield are given in figure 12. Using this expression for secondary yield, Purvis, et al. 
(ref. 11) find 



for - 0 


(A9a) 


where i^ erfc|/kTg/^ is the fifth repeated integral of the error function complement. 
Equation (A9a) depends on in a simple manner. It is more convenient to express 
this current density by 


/ 


where 




for 



(A9b) 


is a constant for each combination of electron environment and material properties. 

The appropriate expression for cp > 0 is quite complex. However, for >0, 

S o 

the emitted secondaries will be attracted back to the surface so that their total contribu- 
tion to the current balance will become small as (p becomes positive. For simplicity, 

O 

the expression chosen for I in this case is 

O 
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for > 0 
s 


(A9c) 


= i- 


kT„ 


where kT is the temperature of the secondaries. This temperature is expected to be 

a few electron volts. The present analysis uses kT„ = 2 eV. 

s 


Magnetic Field Effects 

Two effects of Jupiter's magnetic field on current collection were considered: 
gyration around the magnetic field and induced V x B potential gradients. Neither is 
thought to be of great significance for probe charging, since they probably will affect 
potential by only a few tens of volts. Whipple (ref. 8) discusses the nature of the effects 
in more detail. 

Gyration around magnetic field. - Charged particles in a magnetic field will gyrate 
about the field lines with a radius of gyration given by 


* 

P 


/2MkT 
y eB 


(AlO) 


where B is the magnetic induction. If p* becomes smaller than the radius of the 

0 

probe, the electron current to the probe will be reduced. Since ions have much larger 
radii of gyration than electrons, their collection is unaffected. Present calciilations in- 
dicate that at about 3 Jovian radii the thermal electron current will be reduced to 90 per- 
cent of its value with no magnetic field; and at about 1. 4 Jovian radii, to about 55 per- 
cent. Estimates of the effect of this decrease on probe potential indicate that cp^ will 
become less than 10 volts more positive at 1.4 Jovian radii, an insignificant amount for 
charging considerations. 

Induced V x B potential gradients . - The motion of the probe across the magnetic 
field of Jupiter gives rise to a potential gradient across the probe in the V x B direc- 
tion, where V is the velocity of the probe relative to the magnetic field B. Since B 
corotates with Jupiter, this V is the same as the velocity of the probe relative to the 
the corotating plasma (eq. (A7)). For the geometry of figure 3, if the magnetic dipole 
were alined with the Jovian axis of rotation, the gradient would be in the radial direction 
and woxild result in an estimated 10 -percent reduction in the electron current at about 
1. 5 Jovian radii. The potential gradient across the probe at 1 Jovian radius is estimated 
as about 20 volts. Because this is a small effect from the charging standpoint, the 
probe's spin would reduce the gradient below the calculated values. The effect could be- 
come significant for a larger spacecraft, since the 20 volts calculated here represents 
the development of a gradient of about 20 volts per .meter. 

16 



The thermal and photoelectron currents to the forebody and afterbody exterior sur- 
faces are calculated from 

^PHO "^e + + ^s ^PHO 

Expressions for the four righthand terms are given in table I for the conditions used in 
this analysis. The term I. is actually I. or depending on whether the probe is 
effectively in motion (region 1, use or at rest (region 2, use I. 2 ) with respect to the 
plasma. Entries subscripted a refer to afterbody parameters; those subscripted f 
refer to forebody parameters. Some currents in equation (All) are zero for certain 
cases. For example, IpjjQ = 0 on the shadowed half of the probe (the forebody) and 
= 0 on the afterbody side in region 1. Figure 5 depicts the nonzero plasma and photo - 
electron currents to the probe in regions 1 and 2. 
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APPENDIX B 


CURRENTS DUE TO HIGH-ENERGY PARTICLES 

Since some high-energy particles can penetrate into the interior of the probe, charg- 
ing of both the exterior and interior surfaces can be expected. The situation is depicted 
in figure 6 and discussed in the main body of this report. 


Currents to Exterior Surfaces 

High -energy particle currents to the exterior surfaces are due to particles with 
energies e < Cj for the forebody surface and e < afterbody surface. 

Electron current to the forebody surface is given by 

where the is electron flux, is the lowest high energy being considered for elec- 
trons, and is the electron energy required to completely penetrate the forebody in- 
sulating layers (e^ for electrons). The ctmrent due to incident high-energy ions (H'*') is 
similarly defined by 


^Hi ^ “ -^f® 

where the T- is ion flux, €q 2 is the lowest high energy for ions, and is the ion 
energy required to completely penetrate the forebody insulating layers (e^ for ions). 
C\irrents to the afterbody exterior surface that are due to incident high-energy electrons 
and ions are defined in a similar fashion, with A^ replaced by A^ and and 
replaced by and ®al ®a2 respectively, the afterbody ener- 

gies for electrons and ions. 

Emission processes such as electron reflection and secondary -electron emission by 
electron impact are insignificant for these high-energy particles. However, electron 
backscattering and secondary -electron emission by proton impact can be significant. Be- 
cause backscattered electrons are taken to be a fixed fraction of incident electrons, 

%S = 
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where ^ is the backscatter coefficient. Backscattered electrons are assumed to 
escape. 

Secondary -electron emission by ion impact Igp is similarly taken to represent a 
fixed fraction of the incident ion current. However, these secondaries leave the svirface 
with low energies, comparable to the energies of secondary electrons released by elec- 
tron impact (i. e. , a few volts). Thus, 


^SP ~ Vhi 



for (p„> 0 
s 


where 




is the coefficient for secondary -electron emission by ion impact. 


(B4b) 


Currents to Internal Surfaces 


Particle currents to the aluminum structure are, as discussed, caused by particles 
impingii^ from the forward direction with energies s e < ^ and to those imping- 


ing from the aft direction with energies :S c < e_. 

a, 6 


If we maintain the convention of an 
added subscript 1 for electron energy, and an added subscript 2 for ion energy, these 
current contributions can be expressed as follows: 


For electrons from the forward direction, 

%e, f " “ ^f® [^e^®fl^ ■ ^e^®f, Al, p] 

For electrons from the aft direction, 

%e.iU=jV[ra<‘al>-reM 

where there is an additional factor of 1/2 in equation (B6) because these electrons are 
considered to be collected by the flat part of the aluminum structure (area of yrr^ Un- 

/ ?\ V P / 

stead of the hemispherical part farea of 2?rrpj. 
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For ions from the forward direction, 


^Hi, f ■ ~ ^f® [^i^^f2^ ■ ^i^^f , Al, 2^] 

For ions from the aft direction, 

^Hi, Al " ^ h^^a2^ ■ ^i^^e2^] 

Particle current to the interior of the afterbody insulating layers is caused by par- 
ticles approaching from the forward direction with energies ai — ^ Thus, for 

electrons 

^He, a " ■ 2 ^a® f e^®f , Al, 

and for ions 

%i, a = " [^i^®f , Al, 2^ " ^e^®e2^] 

Emission of such particles as backscattered and secondary electrons from interior 
surfaces has been ignored for simplicity. In fact, any such emitted particles are ex- 
pected either to return to their surface of origin or to strike some other interior surface. 
The net effect should be to reduce interior potential differences to some extent. 

It is now possible to write expressions for the high-energy c\irrents identified in fig- 
ure 4 and required for the overall current balance. For the forebody surfaces, the ex- 
terior and interior high-energy currents are 

%e,f = "%e^ ^Hi + 

= j V {[ri(‘02> - r'i('t2>]» + W,f) - [re('oi> ' re^i>](i - ?w} (bid 


20 



^HEP, a “ f ' ^He, A1 ' %e, Al^ 


where 


= J V { [^l<'f2' - Al, 2>] - [re<'fl> - AI, 1>]'^ 

+ ^ [ri<‘a2> - rifre2>] ‘ \ [''e(Sl> ' ''e(Sl>]''} (B12' 


and 


W,f- < 


6p for <p^^0 


etpf 

6„exp I- 

P V kTg 


for (p^> 0 


r 


P = < 


1 for |<}7fj 
e<Pf' 


1 + 


'01, 


for <P^>€q^ 


^<Pf 




exp for < -€qj 


01, 


0B13a) 


(Bl4a) 


is the factor which accounts for effects of <p^ on high-energy electron collection for 
large potentials. For the afterbody surface, the exterior and interior high-energy cur- 
rents are 

iHe.a =■ J V {[r’l('02> ' r'i(£a2>]'l * %P.A> ‘ ^<'01) ‘ re(^al>]» * ®'5) 

%EP,f = J V {[’'i('f,Al,2> - ri(<e2>] - ' TeCSl']’') (B‘«> 

where 6jjp ^ and y are the equivalents of 6jjp j and j3; that is, 
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(B13b) 


and 


^HP,a“ 


6p6xp 



for <p„> 0 

a 


1 for|,?.JS£oj 


r = 



for Va>£oi 



for < -£„i 


(B14b) 


These equations are summarized in table n. 



APPENDIX C 


PROBE RESISTANCES AND PENETRATION ENERGIES 
Forebody and Afterbody Resistances 

Both the forebody and afterbody insulators are composed of different insulatir^ lay- 
ers, which are treated as resistors in series. The forebody and afterbody resistances 
Rj and R^^ were calculated from bulk resistivities. The thicknesses of the layers were 
based on information from Ames Research Center personnel and references 6 and 7. 

Bulk resistivity was based on information from Ames Research Center personnel (carbon - 
phenolic heat shield) and from measurements made at the Lewis Research Center (fiber- 
glass honeycomb heat shield). The bulk resistivity of Mylar was assumed to be the same 

1 7 

as that of Teflon and Kapton, about 10 J2-cm. The insulated area was taken to be the 
same as the hemispherical surface area in all cases. The calculation is illustrated in 
figure 13. 


Penetration Energies 


The penetration energies are based on Boeing calculations of the probe's structural 
shielding (fig. 14) and on the range and energy loss tables of Berger and Seltzer (ref. 12) 
(for electrons) and Barkas and Berger (ref. 13) (for ions). For the present calcula- 
tions, the items in figure 14 were grouped into forebody insulators, aluminum structure 
(including the urethane foam), and afterbody insulators; and penetration energies were 
calculated by using the range. and energy loss tables. Table HI summarizes the calcula- 


tions. For the ions, the calculated energies ^ 2 ^e2 greater than 100 MeV 

(table in). Because the environmental model cannot predict fluxes for energies above 


100 MeV, ion fluxes for ^ 2 “ MeV were used, and it was assumed that no ions 
escape the probe completely. 
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TABLE I. - THERMAL PLASMA AND PHOTOELECTRON CURRENTS TO PROBE EXTERIOR SURFACES 





































TABLE m. - PROBE DEPENDENT PARAMETERS USED 


Parameter 

Symbol 

Value 

Basis 

Forebody surface area 

Af 

1.27 m^ 

rp = 0. 45 m (refs. 6,7) 

Afterbody surface area 

\ 

1.27 

r = 0.45 m 

Maximum secondary -electron yield 


4.8 

Ref. 11 


(electron impact) 





Energy to attain 


175 eV 

Ref. 11 


Backscattering coefficient 


0.05 

Estimated 

(ref. 11) 

Secondary emission-electron yield 

«p 

1 

Estimated 

(ref. 8) 

coefficient (ion impact) 





Temperature of emitted secondary 

kTg 

2 eV 

Ref. 8 


electrons 





Temperature of emitted photoelectrons 


2 eV 

Ref. 9 


Resistance of forebody instilators 


4.09x10 D 

Estimated 

(appendix C) 

Resistance of afterbody insulators 

Ra 

3.25xl0^^D 

Estimated 

[appendix C) 

Lowest high energ^y: 





For electrons 

^01 

0.06 MeV 



For ions 

^02 

1 MeV 



Energy to penetrate forebody insulators: 





For electrons 

^fl 

15. 1 MeV 

Estimated (appendix C) 

For ions 

^f2 

99. 5 MeV 



Energy to penetrate forebody insulators 





and aluminum structure: 





For electrons 

^f,Al,l 

16. 5 MeV 



For ions 

^f,Al,2 

^123. 6 MeV 



Energy to penetrate afterbody insulators: 





For electrons 

^al 

0. 7 MeV 



For ions 

^a2 

24.4 MeV 



Energy to pass through probe: 





For electrons 

«el 

17.2 MeV 



For ions 

^2 

^148 MeV 




^100 MeV; model is good only to 100 MeV; Cj ^ 2 = 100 MeV and > 100 MeV used 
(e flux). ’ ’ 




integral particle flux, particles/cm^sec 



5 10 15 20 25 

Distance to Jupiter, Jovian radii 

Figure 1. - Jovian particle fluxes along probe trajectory. 
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Figure 4. - Definition of model currents. 
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5 10 , . 15 20 25 

Distance to Jupiter. Jovian radii 

Figure 8. - Potentials attained by probe afterbody. 



Distance to Jupiter, Jovian radii 
Figure 9. - Differential potentials. 


32 



360 




Rgure 11. - Corotalion effect on ion collection. (Ions are assumed to be collected on only the "front" 
surface of the probe when lv| is greater.than the ion random thermal velocity. > 
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Material 

Forebody 

Afterbody 

Resistivity, 

P. 

Q-cm 

Thickness. 

cm 

Resistance, 

R. 

Q 

Resistivity, 

P. 

n-cm 

Thickness, 

1. 

cm 

Resistance, 

R. 

Q 

Mylar 

loi^ 

0.029 

Z 3x10^1 

10^^ 

0.029 . 

2.3x10^1 

Carbon phenolic 

2x105 

9.33 

147 

— - 



Fiberglass 

10^ 

ZZ 

1.79x10^1 

10^5 

1.2 

9.5x10^0 



1 

4.09x10^^ 



3.25xlrfl 


Figure 13. - Calculation of forebody and afterbody resistances. 
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External 

insulation-> 



Section 

Shielding 

Material 

Density, 

g/cm‘ 

Forebody 

Exterior insulation 
Forebody heat shield 
Fleat shield backup 
Structure 

Multilayer Mylar 
Carbon phenolic 
Fiberglass honeycomb 
Aluminum 

0.11 

7.75 

.05 

.42 

8. 33 

Equipment cover 

Interior insulation 
Structure 

Urethane foam 
Aluminum 

0.16 

.22 

0.38 

Afterbody 

Exterior insulation 
Heat shield 

Multilayer Mylar 
Fiberglass honeycomb 

0.11 

.49 

0.60 


Figure 14. - Probe structural shielding. 
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